Kondo Effect in a Metal with Correlated Conduction Electrons 
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The large-degeneracy expansion for dilute magnetic alloys is extended to account for conduction 
electrons interactions. Particular attention is paid to the renormalization of the hybridization vertex 
which affects the low-energy excitations. As a first example, we calculate the enhanced characteristic 
energy fcaTb in the limit of weakly correlated conduction electrons. The metallic regime with strongly 
correlated electrons is discussed. 
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Correlations among the conduction electrons may 
strongly affect the ground state and the low-energy exci- 
tations of metals with magnetic ions. A typical example 
are doped cuprates with rare earth ions Nd2- a; Ce 3 ;Cu04 
where a novel type of heavy fermion-like behavior was re- 
cently discovered JjJ . Although the thermodynamic prop- 
erties in the normal state closely parallel the behavior of 
standard heavy fermion compounds |2|,|3| the observed 
characteristic energy of the low-energy excitations can- 
not be related to a Kondo temperature estimated from 
the known coupling between the spins of the rare earth 
ions and the conduction electrons 

In this paper, we study the influence of conduction 
electron correlations on the low-energy excitations of di- 
lute magnetic alloys. The latter are characterized by an 
energy scale fc^To much smaller than the characteristic 
energies of the conduction electrons such as the band- 
width D or the local Coulomb repulsion U. We thereby 
tacitly assume that the conduction electron correlations 
do not introduce low-energy anomalies i.e. that the con- 
duction electron properties are smooth on the scale set 
by ksT . This separation of energy scales forms the basis 
of our theory. It allows us to reduce the general problem 
of a magnetic impurity in a metal with strongly corre- 
lated electrons to a form which can be solved in close 
analogy to the well-studied case of uncorrelated conduc- 
tion electrons. Of particular interest is the characteristic 
energy fc^To which in the case of uncorrelated conduc- 
tion electrons depends exponentially upon the inverse of 
the coupling between localized and conduction electrons 
^-^|. This exponential many-body scale is a direct con- 
sequence of the existence of Fermi liquid low frequency 
excitations in a normal metal. Theoretical studies of an 
impurity spin embedded in a one-dimensional Luttinger 
liquid pO|-|r^| , on the other hand, predict a power law for 
the variation of the characteristic temperature with the 
coupling in the limit of large on-site Coulomb interaction 
U. 

The modifications introduced by the conduction elec- 



tron interactions into the low-energy excitations arise 
from the subtle interplay of three different types of in- 
fluences. First, the density of conduction states at the 
Fermi level is changed. Second, the virtual transitions 
between f- and conduction states are reduced. Third, 
the effective spin coupling between the conduction and 
f-electrons is enhanced by the increased number of un- 
compensated spins in the correlated conduction electron 
system. Considering these facts, it is not surprising that 
model studies accounting only for selected aspects arrive 
at rather controversial conclusions concerning the Kondo 
effect in metals with correlated electrons p^l7[ . 

For a microscopic description, we consider an Anderson 
impurity coupled to interacting electrons. The latter are 
described by a one-band Hubbard model. The resulting 
Hamiltonian reads 
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The operators cl (cfc CT ) create (annihilate) conduction 
electrons with momentum k, band energy e& and spin 
(7. The local Coulomb repulsion between two conduction 
electrons at the same site is U. The /„(/m) are the cre- 
ation (annihilation) operators for /-electrons on the im- 
purity site. They are characterized by the total angular 
momentum J and a quantum number m which denotes 
the different states m = 1, . . . , Nf within the iV/-fold de- 
generate ground state multiplet with orbital energy e/. 
All energies are measured relative to the Fermi level. The 
Coulomb repulsion Uf between two /-electrons at the im- 
purity site is assumed to be much larger than the other 
energy scales and therefore we may let Uf — > oo. We 
consider here an orbitally non-degenerate Anderson im- 
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purity putting in Eq. (|l]) m — a and the hybridization 

Coupling Vkma = V. 

The excitation spectra are calculated from the Green's 
functions of the empty state |0 > (i.e. the f° or / 14 con- 
figuration) and the occupied / states \a >, respectively, 
denoted by Gq(z) and G IJ {z) 



where e a = 0, e/ for a = 0, a respectively . The inter- 
actions among the conduction electrons affect the self- 
energies E and E CT which are coupled. They are deter- 
mined by the self-consistent large degeneracy expansion 
in close analogy to the Non Crossing Approximation |lg| ] 
for non-interacting conduction electrons. 

The derivation of the generalized NCA equations pro- 
cedes in two steps starting from a conventional perturba- 
tion expansion for the self-energies E and E CT in terms 
of the bare conduction electron propagators, the bare 
Green's functions for the relevant f-configurations, the 
Coulomb interaction U and the hybridization V . By re- 
ordering and partial summation the series is converted 
into an expansion in terms of bare Green's functions 
for the f-configurations, Coulomb-renormalized propaga- 
tors for the conduction electrons and effective hybridiza- 
tion vertices which account for the correlations among 
the conduction electrons. The resulting diagrams are 
classified with respect to their order in the inverse de- 
generacy in close analogy to the case U = 0. Self- 
consistent summation of the leading terms yields the 
self-energies displayed in Figures [l]a,b for the empty f- 
state and occupied f-states, respectively, while Figure |l|c 
illustrates the general structure of the effective hybridiza- 
tion vertex for the case of two-particle correlations. The 
correlation-induced vertex corrections generally contain 
an n-electron Green's function where In— 1 external lines 
are connected by n Green's functions for empty and oc- 
cupied f-configurations, respectively, via 2n — 1 bare hy- 
bridization vertices. On the relevant low-energy scale 
the variation with energy of the effective vertices and 
of the selfenergies is determined by the structure of the 
diagrams and the Green's functions Go and G m . This 
fact follows directly from the separation of energy scales. 
For the Hubbard model, the dominant contribution orig- 
inates from two-particle correlations. We shall therefore 
restrict ourselves to the self-consistency equations dis- 
played in Figure [l]. 

Neglecting the vertex corrections from the Coulomb 
interaction yields self-energies 

k<7 

E - 0) H = ^E/ #M-04t&OGo(w-0 (3) 
k J 

in close analogy to the case of non-interacting electrons 
§§. In Eq. (§), n F(0 is the Fermi function while 



A a (k, £) denotes the spectral function of interacting elec- 
trons which for U = 0, reduces to S (£ — e^). For non- 
interacting conduction electrons, the self-consistent so- 
lution has three characteristic features: The occupied f- 
spectrum shifts to peak at a value Ef the dominant con- 
tribution to the level shift coming from the continuum 
of charge fluctuations. The resonance in the occupied 
f-spectrum acquires a small width. Finally, the empty 
state spectral function exhibits a pronounced structure 
at uiq = Ef — ksT Q which develops with decreasing tem- 
perature and which sets the scale for the low-temperature 
behavior. This feature is the direct manifestation of the 
Kondo effect reflecting the admixture of /°-contributions 
to the ground state and the low-energy excitations. 

It is obvious (l9) that also for interacting conduction 
electrons the dominant effect of hybridization on the 4/ 1 
configurational spectrum is a shift Ef(U) — e/ = ReTi a of 
the resonance energy which, however, is renormalized by 
the Coulomb repulsion U and its influence on the charge 
fluctuations. The interaction-induced shift, however, is 
rather small and will be neglected in the subsequent 
discussion. The central quantity to be studied here is the 
empty state self energy and, in particular, its variation 
with energy close to Ef which can be deduced form rather 
general considerations. The smooth variation with en- 
ergy of J^k A<r{k, £) implies that the basic analytic struc- 
ture of Eq -* (uj) is not altered as compared to the case of 
non-interacting conduction electrons, the characteristic 
feature being a logarithmic variation in the vicinity of 
the f-energy Ef. The prefactor, however, is proportional 
to the interaction-renormalized density of states at the 
Fermi level p(0). The low-energy scale, /cbTo, i.e. the dis- 
tance between the pole in the empty f-state Green's func- 
tion and the 4/ 1 peak depends on the renormalized pa- 
rameters in the usual exponential way. The vertex correc- 
tion displayed in Figure ^c accounts for the two-particle 
cocrrelations in the interacting conduction electron sys- 
tem. The corresponding contribution to the empty state 
self energy is denoted by Eg (uj). We should like to em- 
phasize that the latter is proportional to V 4 Nf (Nf — 1) 
and hence of the same order in the inverse local degen- 
eracy (1/Nf)° as the leading term E °^ (uj). The same 
classification applies to the 4/ 1 configurational self ener- 
gies E^ and e£^ which are both of the order (1/Nf) 1 . 

In this paper we shall evaluate and discuss the contri- 
bution to lowest order in the effective interaction. The 
conclusions we shall arrive at can easily be extended to in- 
termediate values of U by inserting a more sophisticated 
approximation to the appropriate two-particle t-matrix. 
However, the latter cannot be simply expressed in terms 
of Landau parameters since the microscopic expression 
involves conduction electron energies far from the Fermi 
surface. To study the analytic behavior in the energy 
range of interest, we insert the unperturbed conduction 
electron propagator and obtain 
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9k,k',q = G (uj + e.k- efc+g)G CT (w + ek)n F (ek)n F (-ek+q) 
[n F (ek')G a {uj + ek')-n F (ek'+ q )G rT {uj + ek-ek+q+^k'+q)] 



e f ~> e f 



The occupied state self-energy E^ is expressed analogi- 
cally To make things simpler we calculate Eq 1 '' (lo) 
for T = in the local approximation, hence neglecting 
the momentum conservation in Eq. (0), and obtain 
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Here E 0) (w) = ^ In 
mensionless functions of lo/D. For the square DOS and 



and 1^ and il 2 ^ are di- 



(5) 
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half-filling is given by the following equation 

2 Go(loH-v) 
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A similar expression is obtained for the second integral, 
Iq 2 ^ p3] . Indeed on this first iteration stage of the NCA 
equations we already face an integral equation for the 
empty f-state Green's function ( see Eq. ) . A detailed 
analysis shows that in the energy range of interest the 
integrals vary approximately like ^^(lo) ~ Ailn|o; — 

e/|+ B x and^ 2) (w) (A 2 In \w - e/| + B 2 ) 2 where the 

values of the constants Aj > 0, Bi > 0, i = 1, 2 depend on 
the shape of the empty state spectral function. Detailled 
expressions will be given in a forthcoming paper. For 
a first quantitative estimate of the intgrals 7i,2 we use 
for Go (a;) in Eq. ([jj) the simplest possible form which 
accounts for the charge fluctuations only. This can be 
justified in the Kondo-regime where the spectral weight 
of the charge fluctuations is of order unity. In addition, it 
can be shown that the spin fluctuations do not introduce 
new divergences |l9j. The full NCA iteration procedure 
will be discussed elsewhere [Q. The numerical solution 
for the characteristic temperature fe^To = e / — with 
<x>o given by 



(7) 



is displayed Figure g which compares To (U) as function 
of r for given U and e/. The conduction electron repul- 
sion leads to a net enhancement of Tq(U) as compared to 
the non-interacting case. 

To better understand this result we evaluate the vertex 
corrections at the fixed energy lo = ef — T o (0). This "on- 
shell" approximation which was adopted in |1J] leads to 
the renormahzation 
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Both integrals I\ t2 (wo) are negative, so the former renor- 
mahzation leads to an increase of Tp ( U) whereas the lat- 
ter leads to a decrease ( see also JlT] ). Calculations of 
integrals in Eq. (^) show that the r-renormalization wins 
and To (U) increases with U in agreement with The 
"on-shell" results for Tq(U) which are included in Figure 
H for comparison tend to overestimate the renormahza- 
tion. 

The on-shell approximation in Eq.(^J) may be repre- 
sented as T e ff = r[l + Clip (0)] where C is itself a func- 
tion of T and e f . For e f /D = -0.67 and 0.16 < T/D < 
0.3 we obtain 0.72 > C > 0.54 These values has to 
be compared with the C ~ 1 in the Kondo spin model 
of 0. Analogically we obtain e/ = e/ (1 + CiU/\Ef\) 
with 0.02 > d > 0.016. 

For weak electron-electron interaction the increase of 
To may be understood as resulting from the reduced 
probability of finding doubly occupied and empty lattice 
sites in the correlated conduction electron systems. The 
increased number of uncompensated conduction electron 
spins finally leads to the enhancement of the effective hy- 
bridization coupling. We may expect that for sufficiently 
large U and a half-filled conduction band the virtual tran- 
sition from the f-state to the conduction state will cost 
too much energy inhibiting the To-increase and leading 
eventually to the change in the trend. The analogy be- 
tween the Kondo spin model and the Anderson impurity 
model in its local moment regime is not complete in the 
case of correlated conduction electrons ( see also Jl3| ) . In 
the former case To will increase monotonously with the U 
because of the enhancement of the exchange interaction. 
In the latter case the process is two-staged, it involves 
the formation of a local moment and its interaction with 
the conduction electrons. 

The present results are based on the separation of 
energy scales. We therefore anticipate no qualitative 
changes when using more sophisticated approximations 
for the vertices (1, 2; 3, 4), (Figure |l| c) appropri- 
ate for the strong correlation regime. For a quantita- 
tive treatment of this problem there are two visible ap- 
proaches. One is to introduce summation of the infinite 
subseries of the RPA and the ladder type. In this case we 
introduce in self-energies Sq 1 ^ (lo) and (w) dynamical 
Lindhard functions. The other way is to approximate the 
vertex correction by response functions ( dynamic suscep- 
tibilities ) . This may be done with the use of the local 
approximation |23|-|25|. The work in these directions is 
in progress |l^] . In addition, the influence of conduction 
electron interactions on the spectral properties of mag- 
netic impurities and their dependence upon the doping 
are also interesting topics for future investigations. 

In conclusion, the general NCA equations for the inter- 
acting conduction electrons are obtained and it is shown 



that due to the renormalisation of the hybridization inter- 
action the characteristic energy is increased by the weak 
interactions. 
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FIG. 1. Self-consistent f configuration self-energies and 
contributions to the vertex. The solid, dashed and wavy 
lines represent the dressed propagators for conduction elec- 
trons, occupied and empty f states. The open circle denotes 
the bare hybridization V while open and filled squares are 
the bare on-site Coulomb repulsion and the two-particle ver- 
tex r^2(l> 2; 3, 4), respectively, (a) Empty state self-energy 
Tio{w m ). (b) Occupied state self-energy S CT (iw n ). (c) Contri- 
bution to the effective hybridization vertex, (d) Lowest order 
correction. 
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FIG. 2. T K {U) as function of F for U/D = 0.5 and 
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